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Abstract We investigate the effects of spin on entanglement arising in Dirac field 
in an expanding spacetime characterized by the Robertson-Walker metric. We 
present a general approach that allows us to treat the case where only charge 
conservation is required, as well as the case where also angular momentum conser¬ 
vation is required. We find that in both situations entanglement originated from 
the vacuum and quantified by subsystem entropy behaves the same and does not 
qualitatively deviates from the spinless case. Differences only arise when particles 
and/or antiparticles are present in the input state. 
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1 Introduction 

In recent years we witnessed a continuous tendency of quantum information the¬ 
ory to permeate relativistic physics (see [1] for a review on seminal works along 
this line and [2,3] for more recent achievements). In particular, much attention has 
been devoted to entanglement characterization [4] in relativistic frameworks start¬ 
ing from special relativity [5,6] and moving on to general relativity [7-10]. When 
dealing with curved spacetime it is customary to require flatness in remote past 
and far future in order to unambiguously define particles and antiparticles (often 
the conformally flat Robertson-Walker metric is employed for such purpose) [11]. 
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Taking this approach it has been possible to learn about certain aspects of entan¬ 
glement in curved spacetime [12-16]. In particular, the long known phenomenon of 
particle-antiparticle creation from vacuum has been studied to the end of charac¬ 
terizing entanglement in Refs. [16,17]. In such a context, although Dirac fields were 
analyzed beside bosonic fields, the role of spin has not yet clarified. A first anal¬ 
ysis including spin was carried out in Ref. [18] by assuming angular momentum 
conservation besides charge conservation. 

Here we start from a more general scenario where only charge conservation is re¬ 
quired and we explicitly find the unitary representation of the involved Bogolyubov 
transformations. With this approach we are also able to survey the scenario where 
both angular momentum and charge conservation are required (treated differently 
in [18]). It results that in both cases entanglement originated from vacuum and 
quantified by subsystem entropy behaves the same and does not qualitatively devi¬ 
ates from the spinless case. This holds true also when considering two spin particles 
and two spin antiparticles in input. Moreover, entanglement is not created for in¬ 
put charge equal to 2, while it is created for input charge equal to 1 but only in 
the presence of spin. Finally, a remarkable difference between the cases of charge 
conservation and charge and momentum conservation only appears for having one 
particle and one anti-particle in the input state. 

The paper is organized as follows. In Section 2 we recall the Dirac spinors in 
expanding spacetime and present the most general Bogolyubov transformations 
connecting remote past and far future (input and output respectively) in Section 
3. In Section 4 we construct from the given Bogolyubov transformations an input 
vacuum resulting to only respect charge conservation and evaluate the entangle¬ 
ment generated starting from it. Then, in Section 5 we argue the Bogolyubov 
transformations leading to an input vacuum that shows angular momentum con¬ 
servation besides charge conservation. At the end also the spinless case is revisited 
(Section 6). Although, for the sake of simplicity, the main body of the paper fo¬ 
cuses on the input vacuum, in Section 7 we draw our conclusion by comparing the 
results obtainable for various (pure and factorable) input states. In Appendix A 
the unitary representation of Bogolyubov transformations is discussed. The cor¬ 
responding unitary operator is worked out in Appendix B. Finally Appendix C 
contains the expressions of subsystem entropies when starting from excited input 
states. 


2 Dirac spinors in expanding spacetime 

We start by considering the Robertson-Walker line element (see e.g. [11]) 

= a^(r)(—-|- dx'^ + dy^ -|- dz^), 

where the dimensionless conformal time r is related to cosmological time t by 
T = f a~^(t)dt. Here a(T) is the scale factor determining the spacetime expansion 
rate. We assume flat spacetime as r —>• ±oo. The covariant generalization of Dirac 
equation for spin-^ field i/j of mass m is given by 


( 7 '^{ 9 m + r^)+m)ip = 0 , 


( 1 ) 
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where the curved gamma matrices 7 ^ are related to flat ones through 7 ^ := a ^ 7 ^ 
and 


-Tp = I ^[>. 70 ], 


are the spin connections. Here and below dot denotes the derivative with respect 
to conformal time r. 

Writing 'ip = dv — M)(p, with M = ma^r), in (1) we get 


dfidv^ — 7°M</3 — M^ip = 0, 


( 2 ) 


being the flat metric as opposed to the actual spacetime metric := a ^g^’^ . 
Next, given the flat spinors and (with d =t, i) satisfying the relations 

7°Wd = -iud, 'y^Vd = ivd, 


we set 


^p{ ) .— )j( ) (T-)M^e**’’^, 

with X, p position and momentum vectors in The functions obey the 
differential equation 


+ M'^ ±iM^ =Q. (4) 

Let us define by the solutions behaving as positive frequency modes with 

respect to conformal time r near the asymptotic past {in) / future {out), i.e. 

Ein/out = + ^in/ouV ^^n/out = ma{T -^ -/ + oo). 


Then we can introduce spinors that behave like positive and negative energy 
spinors, respectively, in the asymptotic regions: 


U (p, d, X, r) := N{'y''dv - M)f^ '> {T)ude^, 
V (p, d, X, 7 ) := N{'y''dv - M)f^^'> 


with the normalization constant N = 
conjugation. 


_ 1 _ 

y/2M{E+M) 


and 


denoting the complex 
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3 Input-Output transformations 

Let us assume {Uin, Vm} and {Uout,Vout} to be two complete sets of mode solutions 
for the Dirac equation (2) which define particles and anti-particles in asymptotic 
regions and have corresponding vacua, |0)i„ and |0) out respectively. Physically, 
|0)i„ is the state with no incoming particles (anti-particles) in remote past and 
|0)out is the state with no outgoing particles (anti-particles) in the far future. When 
T —> ± 00 , the spacetime is flat and the dynamics of the field is that of the free 
field. So we have two natural quantizations of the field, associated with two Fock 
spaces. The Dirac field operator can hence be written as 

V’{x,77) = MpX] [a„(p,d)C/i„(p,d,?7)e*P ’^ ± b|'„(p,d)Vj„(p,d,?7)e“*P ’^j 

d 

= I dp^ |^ao„t(p,d)17out(p,d,?7)e*P’'±6))„j(p,d)V'out(p,d,?7)e“*P'’^j , 

•' d 

( 5 ) 

where , bin and aout , bout are annihilation operators of particles and anti-particles 
in the in and out asymptotic regions respectively. Actually atn and differ from 
flout and bout because they do not correspond to physical particles outside the in 
region. However, it is possible to relate the operators of m-particles to those of 
out-particles by Bogolyubov transformation [11]. 

We can find the Bogolyubov transformations relating the input regions with 
the output in the following way. Properties of the equation (4) give the following 
relation for its solutions 

/ih + B^~\p)f^'^KutiT), 

where p = |p|, and and are the Bogolyubov coefficients. 

For r —>• 00 Eq.(6) reads 

A(+)(p)e“*-®‘’“* ± B(+)(p)e*-®““* 

Inserting (7) in (5) and equating the coefficients of curved-space spinors we arrive 
at 


where 


a„(p,d) = A*{p)aout{Tp,d) - ^(-p)&Lt(-P. ^0. 

d' 

^L(-P- =Y1 Pd'di-p)aout{p, d') + A{p)bl^^{-p, d ), 

d' 



Eout ^out “h ^out 
^in ^in ± Min 


A^' \p) 


( 8 ) 


( 9 ) 
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and 


Pdd'iv) = -i 


){p) 


( 10 ) 


^^{Eout + Mout){Ein + Min) 

Note that /3d(i'(~p) = ~Pdd'{p) • The corresponding inverse transformations are 


float (p, d) = Aip)ain{p, rf) + X] Pdd' (-P)^^„(-P. d'), 

d' 

^Lt(-P>'^) = -X]^d'd(-P)am(p,d') + -4*(p)feL(-p,d) . 


( 11 ) 


Hereafter, without loss of generality, the quantity A{p) will be considered real. The 
expectation value of the out number of particles in the in vacuum (i.e., number of 
created particles) with momentum p and spin projection d is 

n^(p,c?) := ^ |/3d,d'(-p)P, (12) 

d' 

where the adjoint of (8) has been used. Similarly, the expectation value of out 
number of anti-particles in the in vacuum (i.e., number of created anti-particles) 
is 

w“(-p,d) := {Oin\blnt{-P,d)bout{-P,d)\0in) = ^I^d',d(-P)l^- (13) 

d' 

Then, the total number of created particles with momentum p (total number 
density) results by combining (12) and (13) and summing up over spin projections 

«-(p) := XI 

d 

From now on, to simplify the notation, we will omit the argument of Bogolyubov 
coefficients by assuming it to always be p, unless otherwise stated. We are legit¬ 
imate in doing that because of the absence of mode mixing in the Bogolyubov 
transformations (8) and (11). Also, we assume that ladder operators a («!) always 
refer to particles with momentum p, while operators labeled by b (fel) always refer 
to anti-particles with momentum —p, unless otherwise stated. 


4 Charge conservation 


The unitary operator (acting on the Fock space) that represents the transformation 
(8) has been derived in Appendix A. By applying it to the out vacuum we will get 
the in vacuum, that is |0p; 0_p)m = b(\0p', 0-p)out- This notation refers to the fact 
that if an observer in the input region detects some vacuum state, an observer in 
the output region will detects some “evolved” state which depends on the dynamics 
of the universe. It results in (see Appendix B) 


2 / Pk Pk 

|0p; 0_p)i„ —A ^|0p; 0_p)o«t -^\ '\p\\r—p)out -^\ \rp'A—p)out 

Ipil-p/out -^1 4p! 4-p/oui “I -^2 


t'I'P) t'f—p}o«t) ■ 

(15) 
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This is the way a vacuum state at in is seen at out, i.e. no longer as a vacuum. 
Notice that the states at the r.h.s. have zero net charge likewise the one at l.h.s., 
however the angular momentum is non zero at the r.h.s. (differently from the 
l.h.s.) due to the presence of states like | \rp-,i-p)out and | tp;t-p)out- The particle 
anti-particle density operator corresponding to (15) in the out region reads 

Qp,-p — 2"^ 0-p)(0p; 0-pl — 2-^|0p; 0_p)(tp;4—p | 

“ 2-^|0p; 0_p)(4,p;t-p | — 2-^|0p; 0-p)(tp;t-p I 
-2^|0p;0_p)(;p;;_p | 2 |0p; 0-p)(np; t-f-p I 

H I tp j 4—p)(tpi i-p I + 2 "^^2 I tpi 4—p)(4-p j t-p I 
+ 2 -^2 I tpj4—p}(tpjt-p I + 2 -^2 I tpi4'-p}(4'pi4'-p I 
_ 2^*^ %%-^%% ! tp;;_p}(np;n-p I + ^^1 ip;t-p}(;p;t-p 
+ 2 ^^2 14'pjt-p}(tpjt-p I + 2 '*'^2 I 4-pit-p}(4'Pi4'-p I 

_ 2^*^ %%-^%% ! I + ^^1 tp;t-p}(tp;t-p 


,/3n% 




+ 2 —-^ 2 —I tp;t-p)(4-p;4—p I ~ 2/3|.|.--1 tp;t-p}(t4'p;t4—p 

+ ^^1 ip;;-p}(;p;;-p I - ;p;;_p}(np;n-p I 

+ ) I t;p;n-p)(np;n-p i] +h.c.. (le) 

The reduced density operator for particle is obtained by tracing over antiparticle 


(out) _ 


= Tr_, 


(«) 


and it results 

=yl"‘|0p}(0p| -I- + |^-|-tl^)| tp}(tp I + 4 'p)(4^p I 

+ (l/3nl^ + t4-p)(t4^p I +tp}(4-p I 

+ + ^-t-t/^u)! 4-p)(tp I • (17) 

The latter can be simplified by accounting for the fact that 

(/Iff^tt + P\iPu) = = 1* > (1^) 

as shown in Appendix A. In the case of only charge conservation, the Bogolyubov 
coefficients are related to the total number of particles n (see Eqs. (12), (13) and 
(14)) in the following way 


A'‘ = K^, |%|" + |%|’ = 5 . 


4 . 4- m 

To evaluate the amount of particle-antiparticle entanglement of (16) we can 


use the subsystem entropy [4] (because the state is pure) 


( 20 ) 
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That leads to (together with relations (18) and (19)) 

s (,r‘>) =-2 (^) log. (^) - 2 (5) log. (s). ( 21 ) 

It is worth remarking that (21) is a concave fnnction of n taking minimum value 
(zero) for n = 0 and n = 4 and maximum value (two) for n = 2. 


5 Charge and angular momentum conservation 


To have angular momentum conservation the coefficients of | \.p\\.-p)out and | tp 
;t-p)oui in (15) would have be zero, however, based on (10), dtt and are non 
zero. This means that we have to remove them from (8) which simplifies to (see 
also [18]) 

a^n{d) = A* aout[d) + 

= -P-d,d aout {-d)+A (d). 


The unitary operator acting on the Fock space and representing such transforma¬ 
tion in derived in the Appendix A. Applying such unitary on the out vacuum gives 
(see Appendix B) 


„/ p;, pr^ pt,pi^ 

I Op; 0_p)i„ — A MOp; Q—p)out -^ \ dpi \—p)out -^ \ 4-1 t—p) oui H | ddpi d4—p) out 

(23) 

Again this is the way a vacuum state at in is seen at out, i.e. no longer as a 
vacuum. Notice however that this time the net angular momentum of the state at 
the r.h.s. is zero likewise that of the state at the l.h.s. The particle anti-particle 
density operator corresponding to (23) in the out region reads 

4?-p = i[-4^|0p;0_p)(0p;0_p| -2A3%|0p;0_p)(tp;d-p | 

— 2A ,flj_-|-|0p; 0_p)(4,p; d—p I + 2A (i-f4,/34,-|-|Op; 0_p}(ddpj dd—p I 

+ dp;4—p)(dp;d-p I + 2^^/3|j_/34,-|.| dp;d-p}(dp;d-p I 

— 2^|/i.^j_| p^^^\ dp;d-p)(ddpid4—p I + A I dp;d-p}(dp;d-p I 

— 2^%|%|"| dp;d-p>(ddp;dd-p I +1%|"|%|"| ddp;dd-p)(ddp;dd-p l] +h.c 

(24) 


The reduced particle state results 

= A"|0p)(0p|+A^|%|^| dp)(dp I+-4^I%|"| dp)(dp l + l%lVttl'l tdp)(ddp |. 

(25) 

Notice that from (12) we have 


1 %!^ = »^^(d), 1 %!"’ = nP^t)- 


(26) 


Then, referring to (14) and assuming n^(d) = n^(d) = '^-'^(d) = ’T-°(d) = '1-/4, we 
obtain the coefficients appearing in (25) solely depending on n, that is 


= I%|'= !%!' = 7- 


4 


(27) 
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To evaluate the amount of particle-antiparticle entanglement of (24) we use again 
the subsystem entropy (20) obtaining 

s (,<-)) = -2 log, ( 1 ^) - 2 (n) log, ( 2 ) , (28) 

which is exactly the same expression of (21), although the states (16) and (24) are 
not the same. 


6 Entanglement in 1 -|- 1 spacetime 


In the case of 1 -|- 1 spacetime all coefficients of one particle and one antiparticle 
states in (15) would be the same due to the lack of dependence on spin indexes. 
Moreover the coefficients of two particle and two antiparticle states would be zero. 
This means that (8) simplifies to 

fflin = -4 Oout + P 

= —PO'Out + A 

The corresponding unitary operator is derived in the Appendix A. Once it is 
applied to the out vacuum it yields (see Appendix B) 

|0p; 0—p)in = A 


|0p; 0—p)out ^—p)out 


(30) 


The density operator in the out region corresponding to (30) reads 

Qp,-P = I 0-p)(0p; 0-p\-2pA\0p-, 0_p)(lp; l-p\+\P\^\lp-, l-p)(lp; l-p|] +h.c., 

(31) 

and the reduced particle state is 

4°“‘)=^^|0p)(0p|+ 1/31^1 Ip) (lp|. (32) 

Referring to again to (14) we can now assume nP = rA = n/2 and Eq.(27) becomes 

^" = ^1 1 ^ 1 ' = ?' ( 33 ) 

Then, the amount of entanglement of (31), thanks to (20) and (33), results 

s(gr')^-(l-|)log,(l-2)-2log,2 (34, 

Also (34) results a concave function of n taking however minimum value (zero) for 
n = 0 and n = 2 and maximum value (one) for n = 1. 

Comparing (21) and (28) with (34) we end up with the following relation 

^spinipP — spinlessip 1 0 ^ 72 ^ 4. 
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7 Conclusion 

In conclusion, we have investigated the role of spin on particle-antiparticle en¬ 
tanglement arising from the vacuum of the Dirac field in an expanding spacetime 
characterized by the Robertson-Walker metric. We have considered two approaches 
in dehning the input vacuum in terms of output states: one requiring only charge 
conservation, the other where both charge and angular momentum conservation 
was required. We have found that, although the resulting output states are differ¬ 
ent in the two approaches, their amount of entanglement is the same and simply 
related to that of spinless case, namely Sgpinin) = '^Sspiniessip-/'^)-! for 0 < n < 4 
the density of particle creation. We can then conclude that, in this framework, 
spin does not play a relevant role in entanglement generation. More generally this 
is not true. When considering input states different from the vacuum it is pos¬ 
sible to hnd qualitative and quantitative differences for the sub-system entropy. 
Its expression for various input states is reported in Appendix C. Upon making a 
comparison, we can summarize here the following results: 

i) For charge equal to 2, i.e. two spin particles or two spin antiparticles, it is 


Sspin (ri) — 0. 

ii) For charge equal to 1 it is Sspin{n) = - (^) log 2 (^) -(f) fog 2 (f )> while 


^spinless{.'^l‘^) d- 


hi) For zero charge with two spin particles and two spin antiparticles, the situation 

is identical to that of vacuum, where Sspinin-) = ‘^Sspiniessi’’^/"^)- 
iv) For zero charge with one spin particle and one spin antiparticle we have: 

a) in case of only charge conservation Sspin becomes function of an additional 
parameter A, besides n, in a way depending on whether the spin are parallel 
or anti-parallel; 

b) in case of charge and momentum conservation Sspin only depends on n and 
it is Sspinin) = 2Sspiniess{iT-) for anti-parallel spin, while S'spm{ri) = 0 for 
parallel spin. 

Hence, provided that the input state is not the vacuum, by measuring entanglement 
we could distinguish between different scenarios. 

Finally, spin may play some role on entanglement degradation once entangle¬ 
ment between modes of different momentum is considered (as pointed out in [18]). 
Hence the irrelevance of spin in entanglement generation starting from vacuum 
should be ascribed to the absence of mode mixing for the considered model of 
spacetime dynamics. As consequence we can argue a relevance of spin, for vacuum 
originated entanglement, in anisotropic spacetime where mode mixing appears. 
This is an avenue to take for further studies. 


A Unitary representation of Bogolyubov transformations 

Below we will derive the unitary operator (acting on the Fock space) that represents the 
transformation (8) (and consequently of (22), (29)). In fact the coefficients of this latter, once 
arranged in matrix form, can be regarded as an element of a matrices group [20]. 

Let us start considering fermionic mode operators fi, satisfying the algebra 



{fi,fj} = 0. 


(35) 
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For our purposes fi := aoutit), h '■= Uoiit(4), h ■= fcout(t)i and U ■= 6out(4-)- Let us further 
define the unitary operator (fermionic ‘squeezing’ operator) 


U := exp 




(36) 


by means of complex coefficients that have to be intended as entries of an antisymmetric 

matrix 6 with null diagonal elements. 

We want to obtain the input (m) operators by the action of lA on the output {out) ones. 
For example we want ain{'\) = lA aout{'\)lA^• Let us then compute 


where 


U^j — fj + — [LuJj] + + • • ■ 


We will have 


2 ^ 
*1^2 


i 

[Lu[Lu[Lu,^j]]] = Y. 
[Lu[Lu[Lu[Lu,^jm= E 


(37) 


To show how to derive such relations we provide the explicit calculation for the first of them, 
namely 


iLuJj] = i E {0ni2VlA’fi] + 

ili2 

— 2 E Vi2’ 

ili2 

= ^ E ®n*2 (f/i {f|.. fj} - K ’ 

*1*2 

= 2 E ®“1®2 ~ 

ili2 

= ^ (E®n4“Efl2®7 

\ n *2 

i 

The relations (37) give 

lAtjlA^ =fj + — ^ ^^ 2*^1312 ^*: 

i ill ii\i2 **l*2*3 

+ ^ X!. ^*2®*^2313 ^*4*3 ^1514 ^*5 + ’ ’ ' 


5! .. ^ 

*2 *3 ^4 
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Swapping the indexes on 6 and renaming some of them allows us to write 

UfjU^ — ~ 2\ ^3-i-l^iix ^ ^*2^1 ^^2 i ^ ^*2 

i ’ ’ *1 ’ ni2 ’ il«2«3 

~ ^ 55 ^i 2 n ^^2*3 ^ 14*3 ^*4*5 ^zis ^ 


il«2«3M«5 


= ^ {Hi fi + f/) ■ 


(38) 


To determine the coefficients fiji and i/ji we consider the polar decomposition 0 = O\0\ where 
© is a unitary matrix and \6\ is defined through the relation 6^6 = |0p. Hence, we have 

n 

E = E(l®l)k(l®l)?2i = = m^)ii . 

ni2«3 *2 ^2 

(39) 

and so on. Additionally, we can write 

9ji = = E WiMki^du . (40) 

I I kl 

Then, using (39) and (40) into (38), we arrive at 
fiji = [cos(|i9|)]ji , 

Hi = = -[sm(|6»|)|6»ri e]ji . (41) 

kl 


The coefficients fiji and i/j^, have to satisfy the following conditions, coming from the algebra 
of fermionic modes (35), 




— (5j. 


{fi.fi} = 


= 0 


y { {Hlfi'ji + HHji) — ^ij 
I 

+HihO = 0- 


(42) 


Now, we write the Bogolyubov transformations (8) as 


(^•in (T) 


-.4 0 0 O' 


ffil 


(-I-) 


0 ,4 0 0 


f2 


5in(t) 


0 0^0 


f3 


L&in(i)J 


. 0 0 0 .4. 


Lf4j 



0 
0 


0 

-ni 


0 


tt 


0 





f 

. 

LU J 


therefore, the matrices /r and v have to be 


M = 


-.4 0 0 0 ■ 


ro 0 -/3|^-/3|^l 

0 0 0 


0 0 

0 0 .4 0 

5 ^ — 


. 0 0 0 yl. 




(43) 


In passing, it is worth noticing that the block matrix 


p u 
V* p* 


can be regarded as an element 


of a matrices group [20], whose representation in the Fock space is given by hi of (36). 
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The conditions (42), given (43), can be translated into 


and into 


l-4|2 

+ 

\0n? 

+ 

= 1 

|yl|2 

+ 

M 

+ I/3-1.-1.P 

= 1 


+ 

\Pn? 

+ l/^tlP 

= 1 


+ 

I/ 344 P 

+ l/^itP 

= 1 



/ 


■ 



1 



(44) we have 




n!’(t) = nP(;) 

= n“(t) = 


f 

li%i 


— 1 / 54 . 4.1 
= i%i 


Furthermore, from (45) we can derive the following relation 

phPu - niPh = 

P4,t P-l-t 




(44) 


(45) 


/^4t 


= {\0n? + \Pn?V^- 

Coming back to the unitary operator U of (36), we can argue for 9 the following form 


9 = 


0 0 »9l 1?2 

0 0 »94 i?3 

-Si -Si 0 0 

-S 2 -S 3 0 0 J 


(46) 


and impose to be diagonal (because /j. is also diagonal). This leads to 

+ ^1^3 = 0 , + ^1^3 = 0 , 

which give |dip = and |d 2 p = Thus 


= diag 


+ W2?, Vl’^lP + l’^ 2 P 


As a consequence of (41) we have 


/i. = cos 


with 


+ \S2\‘ 


10 0 0 
0 10 0 
0 0 10 
0 0 0 1 


0 = 


Vl’^ip + 


sin + l’^ 2 |- 

TPtF+W 


r 0 0 »92 

0 0 1?4 i^s 

-i?i -Si 0 0 

I-S2 -S3 0 0 


Finally equating (47) to (43) we will get A real and 


0 0 i?i »92 

0 0 1?4 »?3 

-Si -Si 0 0 
-S 2 -S 3 0 0 


(47) 


(48) 


dl = 

arccos A. 


sin (arccos A) 

l?2 = 

arccos A 

^14.’ 

sin (arccos A) 

1?3 = 

arccos A 


sin (arccos ^4) 

^4 = 

arccos A 

/3|f 

sin (arccos ^4) 
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Here, we want to consider Bogolyubov transformations able to preserve angular momen¬ 
tum. They take the following form (see Ref. [18]) 




-.4 0 0 O' 


ffil 


r 0 0 0 -/3|^i 





0 .4 0 0 


f2 


0 0 -/3|^ 0 



&(t)in 


0 0^0 


fa 


0 0 0 


f] 

KiU J 


. 0 0 0 yl. 


[ 14 ] 


0 0 0 J 


[flJ 


therefore, the matrices /i and u have to be 


-A 

0 

0 

0 • 


■ 0 

0 

0 

^n 

0 

.4 

0 

0 


0 

0 

-7t 

0 

0 

0 

A 

0 

5 ^ — 

0 


0 

0 

. 0 

0 

0 

A 


L/^n 

0 

0 

0 


This implies that the matrix (46) simply becomes 

■ 0 0 0 

0 0 0 

^ ~ 0 -'i?4 0 0 

. -i ?2 0 0 0 . 


with |'d 2 | = |'d 4 | and |0| = diag[\'d 2 \]’ As consequence it is 


^ 

W2\ 


and also 


= cos(|i? 2 |) 


ri 0 0 0] 


r ° 

0 

0 


0100 

sin(|i? 2 |) 

0 

0 

1?4 

0 

0010 

, u — - 

|^2| 

0 

-^4 

0 

0 

Lo 0 0 ij 

1-’B2 

0 

0 

0 J 


(49) 


Finally, for the spinless transformation (29), the matrix (48) simply becomes 



0 •& 
-i9 0 


with =- arccoB ^ at unitary operator reads 

sin(arccos JK) ^ r 

U = exp + 1 ?* aoutboutj ■ 


(50) 


(51) 


B Decoupling the unitary operator 14 

In this Appendix we will show how the unitary operator defined in (36) acts on a vacuum 
state. First of all it is useful to recall the following result for the exponentiation of operators 
belonging to the su(2) algebra [21]. 

Given operators K± and satisfying the following commutation relations 
[K+,K_] = - 2 K 3 , [K 3 ,K±] = ±K±, 


it is 


g(7+i^++7--ftr-+73i<^3) = ^r+K+^{\nr^)K3^r_K^ 


(52) 


r± = 


2^ cosh ^ — 73 sinh ^ 




with 
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For our purposes we take 


-^+ “ 2 r ^ 
«li 2 


2r ^ 
« 1»2 


then = r and 73 = 0, with r the eigenvalue of the matrix \6\. It also results /3 = difr, 
Fs = —, r± = tanr. 

Furthermore we take 

^3 = iEKf^2-4/), 

*2 

where I is the identity operator. 

It can be easily shown that with our choice of the operators K± and we are allowed 
to write the exponential (36) as 


I il*2 


exp 

tanr ft ft 

2 p Z^^m2 ii i2 

< (cos^ r) exp 

- ln(cos2 r) ^ r 

> exp 

tan r v—v . 

2 p X^®n*2fn7 


21 12 

1 

L il J 

J 

1112 


(53) 


Hence when applying the unitary operator U to the vacuum we are left with 

tan^ r / 1 


0 — p')out — cos T ^ ^ 


(r)^n! 


^ 0*112^x7 ^"1^2 1 |0p; 0-p)oii.t ■ 


(54) 


Considering Eq.(46) and f| = give, in the most general 

case we considered 


* 1^2 / 

* 1^2 / 

f) 

h^2 


= ’?l«L(t)6Lt(t) +i?2aL,(t)6Lta) +’?4aL(4,)bL,(t) + i?3aL,(i)6Lt(4.), 

2 

= 2(i?2i? 4 - i?ii?3)aL(t)aLta)6Lt(t)6Lta). ( 55 ) 


n 


= 0 , n > 2 . 


The expression (54) can now be readily computed as 

W| 0 p; 0-p)out = cos^ r {1 + 7l«Lt(t)6Lt(t) + ’? 2 aLt(t)feLt (i) + ’^4a7t(4-)^Lt(t) + ’?3aLt(4-)6LtU)] 

k r 

+ —^(’^2194 -i?li93)aLt(t)aLt(4-)6Lt(t)fcLt(-l-)}|0>o«t 

= .4^1 I Op; 0 —p)out + — [7tl Tpi 1"— + 74-1 \—p)out + 7tl \—p)out + 7t I 4'—p)out] 

(56) 

+ ^(7-1 “ 7t74-)l t'i'Pit4-p)out}, 


where we have used the fact that r = + \'& 2 \^ and hence cos r = A. 
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In case of populated in-states we have also to account for the action of the first and second 
term in Eq.(53). Actually the first term can be computed with relation analogous to (55) 
having creation operators replaced by annihilation ones, while the second term (inside curly 
brackets) can be understood by noticing that 


(cos^ r) 


n 


- ln"(cos^ r) 
2 "n! 



-\tp)c 


(57) 


where N is the number of particles-antiparticles in the state Thus, using (53) and the 

position = U\‘)lj)out, we get the following excited states in the out region: 

inp;t4 — p}i7i i 0—p') out + I'P’ ■i—p)out + -Ip', ‘\'—p)out — A/ 54 , 4,1 ^—p)out 

— A0-y-y \ i pi X—p')out + Up;t4 — p')out 1 

I T-tpi T —p)in — ^ 4 ,'!' I 1"pi 0_p)out ~ 4'P i 0—p)out + A| 1"4'P i T —p)outi 

I I'pi 0—p)in = A| i'pi 0—p)out ~ /^4,'I'I T—p)out ~ ^ 4 , 4,1 t''l'pi ■I—p)out5 

|0p; 1'—p)in = A|0p; p)out + /^4,4,l -i-pi TJ'—p)out + ^'I'4,l Tpi 'T'l'—p)out5 (^®) 

I 1''lpi ^—p)in — I T-i-pi 0 —p)out5 

I I'pi T—p)in = A/9'[''['I Op; 0 — p)out — /5'I''I'^'|'4,| Tpi 'I'—p)out + 'I'Pi T—p)out + (A^ + |/5'|'4,P)| tpi T—p)out 

— /5'|''|'/5_|__|^| 4'P i 4—p)out + I 1'4'pi 1"4'—p)out) 

I 4'pi T —p)in — A^ 4 ,'['| 0 p; 0 — p)out + (A + 1 / 54 , 4,1 )l 'I'P’t"—p)out ~ X—p)out + 1"—p)o'U.t 

~ 4'pi4'—p)out ~ A/5-j.j^l 'l'4'pi T4'—p)out- 




In case of angular momentum conservation Eq.(49) and fj = 
^ (4-)} give 


(t), al 


:4): ^LCt 




ili2 

i Ve- ft ft 


= ^2aLt(t)fcLta) + 

2 

= 2^2^4aLt(tWouta)bLt(t)bLta), 


n 


= 0, n > 2. 


(59) 


Then (54) yields 


U\0p;0-p)out = cos^ r {l + [<?2at^t(t)bLt('l') + 

K. p 

5 ^^2-^2't?4a^out(t)^iout('l')^lut(t)^Iut ('I')} |0)out 

= A^l |0p; 0-p)out + — [/5|4,| tpi 4'-p)out + /5|-^| ipi t-p)out] + t4'p; t4'-p)out|, 

where we have used the fact that r = |t52| and hence cosr = A. 

For populated in-states we have to account for the relations (59), their adjoints, and for 
the operator (57) into (53). Thus, using the position |'0)in = U\^l))out■, we get the following 
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excited states in the out region: 


I T'l'pi T'l'— 

I I'pi 't—p)in 
I 'i-pj ^—p)in 

I I'-l'pi T —p)i7T. 
I Tpi 0—p)i7T. 

|0pi ^—p)in 
|0p; 4— p)in 
I I'-l'pi 0 —p)i7T, 


/^1'4'^4'1'I^P’ 'tp'i ■\'—p)out + -^^'I'4'l 'I'pi T—p)oiit + ^ I 1'4'P i p)oiit 

I tp'i't—p)out, 

0_p)oift + ^ I -ipj T—p)ouf ~ "tpi 4—p)oti,t ~ •^/^'t'4,1 't-i'—p) out 

/^Itl I'p) O — p') out i'J'P i T —p) Olif ; 

A\ tp;0_p>ot.t t4'p;t-p>out, (60) 

A\0p'A—p)out H“ ^■\\\ TpiT4—p)out) 

^|0p; \,—p)out — /3j,'|'| 4^p; p)oiif, 

I T'l'pi 0 —p)oiit- 


Finally, in the spinless case by referring to (51) we take 


-.t ui 


K— — — '0 (loutbout 


^3 — —{(^out^out + b]^ut^out I), 


where r = It^I. Then applying (52) we get 


Op; 0_p)oiit — *^ { l^p j ^—p)out I Ip; 1—p)oiit I • 


(61) 


and, using the position \'(p)in = l^\'4’)out, also 


I Ip 5 1 —p}in 

|0p, l_p)iT2, 

I Ipi 0 —p)iTl 


/3|0)oiit H“ Ip 1 _p)oTi.t 5 

I Op, 1—p)ot/.t? 

|lp;0 — p}out’ 


(62) 


C Subsystem entropies for excited states 

In the case of charge conservation the entropy of the reduced density operators can then be 
straightforwardly computed from (58), and results as: 

S ^Tr_p ^ I t'l'pi t'l'-p)oTit(t'l'p; t'l'-p ) = —2 ^^ log2 ^^ — 2 ^ —^ log2 ^ , 
s(Tr_p [w|np;0-p>o^.t{t4.p;0-p|^/tj) =0, 

s(Tr_p [u\tp;0-p}out{U0-p\U^]) = - (^t^) l°g2 - (4) '°S2 Q) ’ 

S (tv-p [w|np;t-p>out(tip;t-p |W^]) = - (- 7 -) i°g 2 (“7“) ■ ( 7 ) '°S 2 (^) - 

S (Tr_p [w|0p;t-p)o^t{0p;t-p |W^]) = - (- 7 -) l°g 2 (“7“) “ {j) '°§2 (j) ’ 
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and 


S ^Tr_p ^ I tp; 'I'_p)out(tp; t—p I ~ log2 


( 1 -A) 


n(4 — n) 


16 


/ , , n(4 — n) 

(1 >>) 4 log2 

/ , , n(4 — n) 

[l + Vl (1 4 

/ , , n(4 — n) 

yi (1 A) ^ iog2 

/ , , n(4 — n) 

4 J 


s(Tr_p [w|4.p;t-p>o-tap;t-p |wt]) =-log2 


rx(4 — n) 


16 


/ n(4 — n) 

/l-A ^ ^ ^ log2 


^n{4-n) 

1 + V 

/l-A ^ ^ ^ log2 

1-V 

1 

1 


where A is a parameter entering into play because the modulus of beta coefficients can be 
directly related to the particles density only via 


|/3nr = (i-A)-, 


l/3nl =^ 4 - 


0 < A < 1 . 


In case of charge and angular momentum conservation, Eq.(60) provide us with 
S ^Tr_p ^ I t-i-p; t 4 —(t'i'P; ti—p I ^ ^ ^ log2 ^ ^ ^ — 2 ’ 

s(Tr_p [w|np;0-p>o„t(t-;p;0-p|w1‘]) =0, 

S (Tr_p [w I tp; 0 _p)o^t(tp; 0 -p|Wl‘j^ = - (“Y") (~ 7 ~) “ (7) '°§2 (^) . 

s(Tr_p [w|np;t-p>o.t(np;t-P l^t]) = - (^) l°g2 (^) - O '°g2 (^). 

S (Tr_p [u |0p;t_p).„t(0p;t_p I W^]) = - (1^) loga (^) " [j) (^) , 

s(Tr_p [w| tp;t-p>out(tp;t-p Ipf]) = 0 , 

5(Tt_p [u\lj,-,t-p)out{lp-A-p 1"^]) =-2(^)log2(^) -2(J)log2(J). 

The main difference with respect to the previous case is that here entropy only depends on 
the particle density n. Furthermore, there is no entanglement creation from one particle and 
one anti-particle with both spin up (or down). This was expected because of the conservation 
of angular momentum. 

Finally, concerning the spineless case, as consequence of (62), we have 
S (Tr_p [w|lp;l_p>o„t(lp;l-p|Wt]) = - (-^) l°g2 “ (f) '°®2 (|) . 

s(Tr_p [w|lp;0_p)o^t(lp;0_p|Wtj^ =0, 

s(Tr_p [w|0p;l_p)o^t(0p;l_p|Wtj^ =0. 
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